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Abstract 



I. INTRODUCTION 



The central quantity in a linear response analysis of a thermal field theory is the po- 
larization tensor and in particular its infrared behavior. It determines the spectrum of the 
quasi-particles and their dispersion laws, and also the dynamical and static screening of 
external fields. In the physics of the very early universe, which is filled with a hot plasma 
of various elementary particles, the gravitational polarization tensor is also of interest. It 
describes the response of the plasma to metric perturbations and its infrared behavior de- 
termines the dynamics of large-scale cosmological perturbations. 

In Ref. [1], one of the present authors has calculated the leading temperature contri- 
butions to the gravitational polarization tensor II"^ 1 ' of a collisionless ultrarelativistic (i.e. 
effectively massless) plasma for temperatures T <C mpi anc k, where perturbation theory be- 
comes applicable. Because the underlying effective action turns out to be conformally invari- 
ant, n a/3Ail/ can be calculated with flat-space momentum techniques and the corresponding 
tensor on a curved background space-time with vanishing Weyl tensor is obtained by a simple 
conformal transformation. In Ref. [2] this has been used to establish self-consistent equa- 
tions for linear perturbations of Robertson- Walker cosmological models and exact analytic 
solutions have been constructed. 

In the ultrarelativistic, collisionless case, the leading temperature contributions turn out 
to have a universal form, only the overall normalization changes in accordance with the 
magnitude of energy density. Some subleading contributions, where it starts to make a 
difference which particles constitute the plasma, have been obtained in Ref. [3]. 

In the following, we are interested in the corrections to the leading temperature terms 
when weak self-interactions are switched on. We take the simplest 4-dimensional model, 
scalar fields with A</> 4 interactions, which allows us to carry the calculations through re- 
summed two-loop order. The emphasis of this paper is on the technical aspects of finite- 
temperature perturbation theory such as the need for resummation in order to extract the 
contributions proportional to A 3 / 2 . It will turn out that in the presence of gravitational 
interactions it is no longer sufficient to resum a local thermal mass term only. As in the 
Braaten-Pisarski resummation scheme for finite-temperature gauge theories, it becomes nec- 
essary to include nonlocal vertex corrections which are generated by "hard thermal loops" 
[4]- 

The implications of our results for self-consistent cosmological perturbations will be dealt 
with at length in a separate, forthcoming paper [5]; some first results have already appeared 
in Ref. [6]. 

II. THE THERMAL GRAVITON POLARIZATION TENSOR 

As thermal matter we take massless scalar particles with quartic self-interactions and 
interactions with the gravitational field according to the Lagrangian 1 



1 0ur conventions are those of Ref. [7]. In particular the metric signature is taken such that 

goo > o. 
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(2.1) 



If r denotes all contributions to the effective action besides the classical Einstein- Hilbert 
action, the energy- momentum tensor is given by the one-point 1PI vertex function 
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and the gravitational polarization tensor by the two-point function 
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From the last equality it is clear that II"^ 1 ' describes the response of the (thermal) matter 
energy-momentum tensor to perturbations in the metric. Equating II"^ 1 ' to the perturba- 
tion of the Einstein tensor gives self-consistent equations for metric perturbations and, in 
particular, cosmological perturbations. 

With massless scalars and £ = +1/6, the matter part is conformally invariant, that is 
T[g] = T[Q 2 g] — apart from the conformal anomaly which like other renormalization issues 
can be neglected in the high-temperature domain on which we shall concentrate. 

This conformal invariance is crucial for our prospected applications for two reasons. 
Firstly, it allows us to have matter in thermal equilibrium despite a space-time dependent 
metric. As long as the latter is conformally flat, ds 2 = <t(t, x)[dr 2 — dx 2 ], the local temper- 
ature on the curved background is determined by the scale factor a. Secondly, the thermal 
correlation functions are simply given by the conformal transforms of their counterparts on 
a flat background, so that ordinary momentum-space techniques can be employed for their 
evaluation. 

The energy-momentum tensor through two-loop order is given by the diagrams of Fig. 1 
which are easily evaluated in flat space-time, e.g. within the imaginary-time formalism by 
choosing periodic boundary conditions for the scalar field in imaginary time with a period 
equal to the inverse temperature [8], 



9=V 



90 ~ 48 + 



T 4 {K 8 ° - Vr) 



(2.4) 



With 7] — > g = arj one has T f 



8°. 



as the normalized velocity field 



of the plasma, so that the temperature in curved space-time is obtained by T — > a 1 / 2 T. 




+ 1/2 ^w^^^^) + 1/4 




FIG. 1. One- and two-loop diagrams contributing to the energy-momentum tensor of A</> 4 the- 
ory. Wavy lines denote external gravitons and straight lines scalar particles. 



3 




+ 1/4 




+ 1/2 



+ 1/4 





+ 1/2 



+ 1/4 





+ 1/4 



+ 1/8 





FIG. 2. The one- and two-loop diagrams of the gravitational polarization tensor. 



The gravitational polarization tensor is an inherently nonlocal object, already in flat 
space-time, as is apparent from its diagrammatic expansion, shown in Fig. 2. On curved, 
but conformally flat, space it is then determined by 



IW/3<£,y) = o-(x)o-(y)Il^ ctl s(x - y) . (2.5) 

g=ar) g=t] 

Denoting by n^afiiQ) the Fourier transform of II MI/a ^(a; — y) , the numerous com- 

9=V 

ponents of the polarization tensor are restricted by diffeomorphism invariance through the 
Ward identity 



and by conformal invariance through the Weyl identity 



»riw(Q) = -T af) . 



(2.6) 



(2.7) 



Using if", iifx = and = (Q ,q), one can build 14 tensors to form a basis for 
Il tJ,l ' al3 (Q) = pYj\=i c i{Q)T- lval3 (Q), see Table I. The above identities, however, reduce the 
number of independent structure functions to three. Choosing (p = T 00 ) 



A{Q) = il 0000 {Q)/p, B(Q) = n / (Q)//>, C(Q) = U^{Q)/p 
the C1...14 are determined by the linear combinations given in Table II. 



(2.8) 



III. ORDINARY PERTURB ATIVE CALCULATION 

A perturbative expansion in Feynman diagrams makes sense for temperatures T <^ 
G~ x l 2 ~ mpianck- In this regime higher-order diagrams generated by additional graviton 
lines are correspondingly suppressed. Apart from this restriction, we shall concentrate on a 
high-temperature expansion in the sense of T ^> Q , \q\(= 5), since the typical length scale in 
the theory of cosmological perturbations is given by the horizon scale ~ (GT 4 )" 1/2 > T _1 . 
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The one-loop diagrams for the energy- momentum tensor and the gravitational polar- 
ization tensor describe collisionless matter. Their leading high-temperature contributions 
have been obtained first in Ref. [1] and the complete underlying effective action has been 
found subsequently in Ref. [9]. Here one may also include a thermal graviton background, 
since the leading-temperature contributions are proportional to the number of spin degrees 
of freedom, so gravitons are as important as any other thermalized matter. The one-loop 
result for the high-temperature gravitational polarization tensor is in fact determined by 
only one independent structure function, since it is constructed from a totally symmetric 
tensor 

K^K^Kfs r = r d 3 k 

Tk~ v hL.J (2tt) 3 ' 



k K\K - Q) 2 



Kq =2izinT ' 



according to 



n (1) 



-yap- 



(3.1) 



Therefore the functions B and C are directly related to the momentum-independent energy 
density, to wit, 

£(1) = _1 C« = 0, 



and all the nonlocalities reside in 

A {1 \Q) 



i 1 5 
u> artanh 

uj 4 



(3.2) 



with uj = Qol q. 

This function has a logarithmic branch cut between uj = ±1. The discontinuity along 
this cut is similar to the one found in the polarization tensor of gauge theories where it is well 
known to correspond to the phenomenon of Landau damping. Indeed, in the application of 
the above results to the theory of cosmological perturbations [2,10] it has been found that the 
discontinuity of A is responsible for the collisionless damping of subhorizon perturbations. 

At two-loop order one begins to see the effects of self-interactions of the thermal matter. 
If A ^ \/~GT ~ T I mpianck-, which we shall assume, only the matter self-interactions are 
important. A straightforward (albeit somewhat tedious) evaluation of the diagrams of Fig. 2 
yields 
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(jO 



uj artanh- 



(jO 



(jO 



or 



(3.3) 
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(3.4) 



C< 2 >(0) 
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3(a> 2 - if (uj artanh- - 2(uj 2 - l)(3u> 2 - 2)uj artanh- + 3a> 4 - 4a> 2 
V uj / uj 



(3.5) 



In addition to a logarithmic branch cut between uj = ±1, whose discontinuity is now 
more complicated, there are also simple poles in A at uj = ±1. 
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(a) 



(b) 



FIG. 3. Two examples of infrared divergent graphs beyond two-loop order. 



IV. RESUMMATION 

If one now attempts to go further, one encounters a problem with ordinary perturbation 
theory starting at three-loop order. Some of the diagrams, e.g. those shown in Fig. 3, are 
infrared divergent. 

The scalar self-energy subdiagrams in Fig. 3a are proportional to XT 2 and their repeated 
insertion generates a chain of massless propagators all with the same momentum. This type 
of divergence and its treatment are well-known [11]. The interactions of the (massless) scalar 
particles with the plasma generate a thermal mass which needs to be resummed. However, 
in addition to repeated insertion of thermal masses there are also 1PI vertex subdiagrams 
which can cause infrared divergences. In Fig. 3b the vertex correction has a similar effect as 
the second mass insertion in Fig. 3a — both make the respective diagram infrared divergent. 

It is therefore necessary to include the (generally nonlocal) vertex corrections when curing 
perturbation theory by a resummation of thermal masses. This can also be understand by 
noting that a simple local mass term would break conformal invariance. But apart from 
the zero-temperature conformal anomaly the full energy-momentum tensor has still to be 
traceless for conformally coupled scalars, which shows that thermally induced masses are 
not equivalent to ordinary rest masses, not even for scalar fields where the thermal mass is 
independent of momentum. 

A systematic way to improve perturbation theory when vertex corrections are of equal 
importance has been worked out first for high-temperature quantum chromodynamics by 
Braaten and Pisarski [4]. It has been found there that it is necessary to resum all the 
subdiagrams which are generated by loop momenta ~ T (termed "hard thermal loops"). In 
A</> 4 theory without gravitational interactions, the only hard thermal loop is given by the 
scalar self-energy diagram, Fig. 4a, which yields a constant mass term m 2 = XT 2 . Clearly, 
when probing momentum scales < y/XT, this has to be resummed. In the presence of 
external gravitational fields, vertex diagrams like the ones in Fig. 4b are of comparable 
magnitude, however they turn out to be nonlocal. 

In Ref. [12] it has been shown how to construct the entire generating functional AS[<f>,g] 
for hard thermal loops in the presence of gravity. In the case of A</> 4 theory this becomes 
particularly simple [13]. The hard thermal loop effective action can be written as an integral 
over a forward scattering amplitude J[K;<f>,g] according to 
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(a) (b) 

FIG. 4. (a) The scalar self-energy which gives rise to the thermal mass, (b) The general form 
of the non-local vertex diagrams which have to be resummed when going beyond 0(A). 



AS[<M = j^y 3 J d*K6(K 2 )6(K )n(K )J[K;<l>,g]. (4.1) 

where n(K ) = (exp(K /T) — l) -1 denotes the Bose- Einstein factor. It then turns out that 
all the dependence of J on the metric can be eliminated by choosing the diffeomorphism 
gauge condition 

K»g, v = K^, (4.2) 

whereupon 

J[K;<f),g] = -6A J d 4 x<f) 2 (x). (4.3) 

Assuming an asymptotically conformally flat space, one can reconstruct the generally co- 
variant function J in terms of the set of null geodesies y^(x, #), y^(x, 0) = x^, whose tangent 
vector y^(x,8) — > for affine parameter 8 — > — oo. The solution is 



(4.4) 



The integral over the geodesic divergence is clearly a highly nonlocal quantity and is what 
reconciles the appearance of a thermal mass for the scalars with invariance under Weyl 
rescalings of the metric. 

The leading high-temperature contributions to the graviton self-energy and higher vertex 
functions with only external graviton lines can also be summarized in a similar (but more 
complicated) generating functional [12]. However, we shall not need it in our resummation 
program since we assumed A ^> (T/mp; anc fe) 2 , by which we can neglect higher-order diagrams 
with internal graviton lines. 

A resummation of all the hard thermal loops summarized by AS[</>, g] is formally achieved 
by writing / d 4 xC = S = S res . - S counter , with S res . = S + AS and S counter = AS. The Feyn- 
man rules now include the thermal mass of the scalars and the additional hard-thermal-loop 
vertices which have two scalar lines and an arbitrary number of graviton lines. Overcount- 
ing is avoided by subtracting them as counterterms, as which they are treated as one-loop 
objects. [Equivalently, one could consider AS as coming from integrating out all the modes 
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with hard momenta > T after which S + AS is the effective action to be taken when 
integrating out the remaining soft modes.] 

Performing a resummed two-loop calculation, we can collect all terms up to but excluding 
order A 2 , which is the formal order of bare three-loop diagrams. Since the resummation 
introduces y/\ through m = \/\T , we can go up to order A 3 / 2 , where the plasmon effect 
enters for the first time. [Because the bare two-loop diagrams are infrared convergent, 
there are no contributions of order A 1 / 2 .] At two-loop order one has to include the thermal 
counterterms. Here it turns out that the diagrammatics is considerably simplified if the 
hard-thermal-loop mass \/\T in AS is replaced by its resummed value, which is given by 
the gap equation 

™ 2 = ~ 12X t V2 1 2 = " - A3/2T2 + • • • • ( 4 - 5 ) 

Tk K z — m z 7r 

Then thermal counterterms can cancel entire subdiagrams, while the sum of the resummed 
two-loop diagrams is unchanged up to order A 2 . 

A. The energy- momentum tensor up to order A 3 / 2 

When switching to the resummed Feynman rules, the one- and two-loop diagrams for 
the energy momentum tensor of Fig. 1 are replaced by the ones shown in Fig. 5a-c, where 
a dot indicates a dressed propagator or a dressed vertex, respectively. In addition there are 
two one-loop diagrams involving thermal counterterms (denoted by crosses), 5d and 5e. The 
mass counterterm in diagram 5e cancels diagram 5b completely (using Eq. (4.5)), bringing 
down the total number of diagrams to two. The vertex counterterm in diagram 5d reduces 
the full vertex in diagram 5a to a bare vertex, so that the resummation of the vertices turns 
out to be superfluous in this case, for there is no hard-thermal-loop vertex which could 
modify the vertex in diagram 5c. 




(d) (e) 




FIG. 5. The resummed energy-momentum tensor. The dots denote the dressed propagators 
and the dressed vertices, whereas the crosses indicate the thermal mass and vertex counterterms. 
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The net result thus reproduces the one obtained by simple ring resummation [8] and 
reads 



g=v 



tk K 2 -m 2 + 6 ^ \hcK*-m 

'I _ A + ^! + (X 2 )} T 4 U 6 y _ \ 
48 12tt J V M ' 

which is traceless as it should be. 



(4.6) 



B. The gravitational polarization tensor up to order A 3 / 2 

In contrast to the energy- momentum tensor, the resummation of the nonlocal hard- 
thermal-loop vertices turns out to be crucial to obtain the complete result for the gravita- 
tional polarization tensor up to and including order A 3 / 2 . A naive ring resummation of a 
local mass term for the scalars would lead to a violation of the Weyl identity (2.7). 

The diagrammatics for the two-point function at resummed two-loop order is given in 
Fig. 6. As before, the two-loop diagrams which correspond to mass insertions can be canceled 
by thermal counterterms when using Eq. (4.5). However, this time the hard-thermal vertices 
do not cancel out completely, except the one in diagram 6b. In diagram 6a, there is only 
a partial cancellation with the vertex counterterms, and in the two-loop diagrams 6e-g the 
dressed vertices are untouched at this order. 

In the application to the theory of cosmological perturbations, we shall need in particular 
the low-momentum limit Qo,q ~ \[GT 2 <C y/\T = m. In this limit the results become 
independent of the coupling constant £ and read explicitly 

Raf3is{Q) = 

= -\f R A* A P {K ■ P Va/} - K a Pp - P a Kp) {K ■ P VlS - K^Pe - P^K 6 ) 
~\^ R ^kK 2 (rj^rjug + Tjasrifr - Tj a/} r] lS ) 

+ \^f R A* Ap {K ■ P Va/} - K a Pp - P a Kp) £ A K A P {K ■ P VlS - K^P 6 - P^K 6 ) 
-^A (rj^-q/ss + VasV/Sj ~ Va/sVjs) (j£ Aifj 
-^Vais Ajf) ^ A K A P {K ■ P VlS - K^Pe - P^K 6 ) 

~\ Xr] ^ (f K Ak ) f K AkAp ( K ' Pr]a P ~ KaP ? ~ PaK ^ 
-^V*f3Vi6 m4 ^ K A*: Ap 

-^Am 2 ^ A^ Ap) £ A K A P (K ■ P VlS - K^P S - P-yK s ) 
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*m%s AkAp^ £ A K A P (K ■ P Vct/s - K a P p - P a K p ) 

18A 2 (jf, AkAp^ A K A P (K ■ P Vct/s - K a Pp - P a K^ 
x£ A K A P (K ■ P VlS - K^Ps - P^Ks) , 

J K 



(4.7) 



where A K = 1/(K 2 - m 2 ) and P = K - Q. 

With this expression one can explicitly verify that the Ward and Weyl identities are 
fulfilled. As a consequence the tensorial structure is as given in Table II, with three inde- 
pendent structure functions. Evaluating the integrals involved as described in the Appendix 
leads to our final result 

, 1 5 

A = u> artanh 

uj 4 



+ 
+ 



5A 
8?r 2 
5A 3 / 2 

8tt 3 



' / 1\ 2 1 L0 2 

2 (jj artanh— — w artanh 

V <J0 / (jO (jO 2 — . 

3 (a; 2 — 1 — uj\/ (jo 2 — l) fa; artanh— J 



+6 uVa> 2 - 1 - a; 2 



(jO 



Vu 2 - \, 



u> artanh- 



(jO 



(jO 



+ 3^ - + 6- 



(U 2 - If I 2 LU 2 -1 



ZloVuj 2 - 1 + 3a; 2 



(4.8) 



B 



+ 



-1 + 



15A 3 / 2 
8?r 3 



/ 1\ 2 1 
-(a; — 1) I u> artanh— ) + (2a; — l)a; artanh u> 

{lo(lo 2 - if 2 - (lo 2 - l) 2 } u> artanh— 



+2 {(a; 2 - l) 2 - a;(a; 2 - 1) 3/2 + loVuj 2 - 1 - a; 2 } u artanh- 



2a;Vu> 2 - 1 + w(a; 2 - 1) 3/2 - a; 4 + 4a; 2 



(4.9) 



C 



+ 



8tt 2 
15A 3 / 2 
16?r 3 



/ i\ 2 i 

3(a; 2 - l) 2 a; artanh- - 2(a; 2 - l)(3a; 2 - 2)a; artanh- + 3a; 4 - 4a; 2 

3(a; 2 — l) 2 (a; 2 — 1 — u>\/u) 2 — l) fa; artanh— J 
-2(a; 2 - 1) {3(a; 2 - l) 2 - 3a;(a; 2 - 1) 3/2 + 3aVa> 2 - 1 - 3a; 2 + l} u artanh- 
-3a;(a; 2 - 1) 5/2 + 6a;(a; 2 - 1) 3/2 + loVuj 2 - 1 + 3a; 6 - 15a; 4 + 16a; 2 ] . (4.10) 
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(e) (f) (g) (h) 




FIG. 6. The resummed gravitational polarization tensor. 

A new feature of the contributions proportional to A 3 / 2 is that in addition to poles at u> = 
±1 and a logarithmic branch cut in between there are now terms involving root singularities 
and a corresponding branch cut. The singularities themselves are in fact spurious, since 
the perturbative expansion breaks down for \u>\ — 1 < A. As found in a somewhat different 
context in Ref. [14], the singularities at u> = ±1 are removed by the thermal masses of the 
scalar particles, although the branch cut for \u>\ < 1 persists. 

Consider the first diagram on the right-hand side of Fig. 6 and its contribution to the 
discontinuity of A = IIoooo/ P-, 

Disc A x = -9(1 - uj 2 ) / dpn'(p) (p 2 - m 2 /2) , (4.11) 

— U)* 

where we have taken the limit Q , q <C m = \/\T but refrained from any further expansion. 
The expansion in A as performed above gives 
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2^ if 5A 15A 3 / 2 \ 
D„cA 1 = 6l(l- W >^l + — -— Jx 

+ ^(l^ + 7r^-(T^) +0(A2) }' (4 - 12) 

which reproduces the discontinuity of A at leading order as well as that of the most singular 
term that ocured at order A 3 / 2 . 

In the limit 1 — o> 2 <C A, on the other hand, one can read off from Eq. (4.11) that 

, x 15u> A 2 / / A \ 

DiscAi -> 6(1 -uj 2 ) — TT rrrexp -J , (4.13) 

V ; 4tt 3 (1 -co 2 ) 2 F ^ V 1 -u 2 ) ' y ' 

so the discontinuous and even singular onset of the imaginary part of A that appeared in 
Eq. (4.12) is illusory. It actually sets in smoothly, but nonanalytically, with all derivatives 
vanishing. Some consequences of this in the context of cosmological perturbations have been 
discussed in Ref. [6]. 

A more complete application of all the above results to the determination of the dynamics 
of cosmological perturbations in the presence of a weakly interacting ultrarelativistic plasma 
will be the subject of a forthcoming publication [5]. 
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APPENDIX: EVALUATION OF THE THERMAL INTEGRALS 

The evaluation of the massive thermal integrals in Eqs. (4.8-4.10) was performed by a 
MATHEMATICA program. The principal steps are the following. 

First a partial fraction decomposition in A; is performed, which gives a sum of terms in 
which A; appears only linearly in the denominators. After a shift A; — > k° + q° where neces- 
sary, the thermal sums are carried out with help of the formula (dropping zero temperature 
contributions) 

T E v^e=±« E ) ( A1 ) 

The sum does not converge as it stands but it occurs only in combinations where the di- 
vergent part cancels and the above formula applies. The energy can take on the values 
+ m 2 or J{k — q) 2 + m 2 according to the loop momentum. In the latter case it 
is convenient to perform a shift of the integration variable k — > k + q such that Ek becomes 
the argument in the distribution function in all summands. 
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The spatial loop-integration is split into an integration over k and the angle 
which is performed first. This leads to the following three types of integrals 



kq/(kq) 




(A2) 



(A3) 



(A4) 



where for the moment we specified to the kinematical region Q 2 = a; 2 (4m 2 - g 2 (w 2 - l))/(w 2 - 
1) < 0. The corresponding formulas in other regions can be obtained by analytic continua- 
tion of the final expressions. 

The high temperature expansion of integrals of the form (A2) is well known [8,15] and 
will not be repeated here. The strategy for the remaining integrals is to subtract off the large 
momentum behaviour of the integrand until the subtracted integrand decreases sufficiently 
for large k such that the Bose- Einstein distribution function can be replaced by T/Ek- For 
the integrals (A3), a suitable subtraction is given by the Taylor-series polynomial in 1/k of 
order I of the logarithm. 

On the other hand, for integrals of typ (A4) the subtractions have been chosen in a 
different way anticipating the subsequent integration. For I = 1 it suffices to subtract the 
^-independent large momentum limit 2artanhu;~ 1 whereas in the case I = 3 an additional 
second order subtraction proportional to 1/-E| was performed. The subtractions themselves 
give rise to integrals of typ (A2) without the energy in the denominator and may be evaluated 
along similar lines as the integrals (A2). The subtracted integral with the replaced Bose- 
Einstein distribution can be evaluated elementarily by introducing Ek as integration variable. 

It still remains to evaluate the order T contribution of the typ (A3) integrals. For this 
purpose we take advantage of the fact that the integrand is even and integrate from — oo 
to oo. Closing the integration contour in the upper half-plane ^m(k) > only picks up a 
cut contribution from the logarithm and a pole contribution at k = im since on the arc the 
subtracted integrand vanishes in the limit \k\ — > oo. 

The final expressions are obtained by replacing the thermal mass by its value given in 
(4.5) and expanding in powers of A 1 / 2 . 
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TABLES 



Tf» v = (vP 7] av + u a rf v ) + U v (vP 7] ali + U a if*) 
Tf^ = u a vP «" u v 

rpapfiv _ a p 
± i — 'I 'I 

Tf» v = U v T) " 3 + u a vP ff v 

Tf^ = U P (Q v TJ afl + T] av ) + Q 13 (ll v TJ afl + ll* T] av ) 

+ U a (Q v if" + rf v ) + Q a {u v if" + ll* rf v ) 

T" p,n/ = Q v u a vP ui 1 + Qi 1 u a vP u v + QP u a ui 1 u v + Q a vP ui 1 u v 
Tf^ = QP Q v Tj " 1 + QP Q» r] ai/ + Q a Q v tjP^ + Q a Q» if" 
Tf» v = Q» Q v u a vP + Q a QP «" v? 
T?t v = (QP u a + Q a vP) (Q v «" + u") 

q^apfiv = Qp Q V u<x + Q a Q V ^ + Q a Qp Q v ^ + Q a Qp ^ 

T^ V = Q a QP Q v 

T^ v = Qp Q v T] a P + Q a QP tf v 

T^ v = {Q v «" + u v ) T] a P + (QP u a + Q a vP) tf v 



TABLE I. A basis of 14 independent tensors Tf^ built from if v ', ui 1 = 6%, and 
<5 M = Q^/q = (w,q/q). 



ci = \Q*A + \Q*B + \C + ±Q* + f& + \ 

c 2 = iq°a + q 4 b + \q*c + £g 6 + i g 4 + ^g 2 -| 

cs = Q 2 { fQ 6 A + iq*b + \Q*c + §§g 6 + H g 4 + f 2 Q 2 - 1} 

c 4 = Q>{lQ>A-±;B-\C + j- 2 Q 2 -%} 

C5 = g 2 { |g 4 A - \q-b - \c + ^g 4 - H g 2 + 
C6 = u { f g 4 A - q>b -\c- ^g 4 - ||g 2 - ^} 

c 7 = U { ^Q°A - IQ*B - \Q^C - §§<? - fg 4 - ^ + |} 

cs = (fg 2 + \) q 2 a + (g 2 + 1) b + \c + ^g 4 + %Q 2 + f- 2 

C9 = (f g 2 + 'i) Q 4 a + (f g 2 + 3) q 2 b + (|g 2 + \) c + ff g 6 + ||g 4 + |g 2 + 1 

cio = (f Q 2 + Q 4 a + (f g 2 + f) q*b + \Q*c + ff g 6 + ff g 4 + f g 2 + \ 

di = « { (-fg 2 - f ) q 2 a - (f g 2 + 1) b - \c - § g 4 - f g 2 - 1} 

C12 = (f g 4 + 5g 2 + 1) a + (f g 2 + 2) b + \c + ff g 4 + ff g 2 + \ 

ci8 = (f g 2 + \) Q 2 a - \q 2 b -\c + £g 4 - l 2 Q 2 - ± 

c 14 = u { fQ*A + \&B + \C - ^g 4 + f A Q 2 - ^} 



TABLE II. The structure of the conformally covariant gravitational polarization tensor 
IL^/p = £!=i CiTr ap in terms of A = iL 0000 /p, B = to^o/p, and C = fi^'/p. 
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